Homeostatic competition drives tumor growth and metastasis nucleation

We propose a mechanism for tumor growth emphasizing the role of homeostatic regulation and tissue stability. We show that competition between surface and bulk effects leads to the existence of a critical size that must be overcome by metastases to reach macroscopic sizes. This property can qualitatively explain the observed size distributions of metastases, while size-independent growth rates cannot account for clinical and experimental data. In addition, it potentially explains the observed preferential growth of metastases on tissue surfaces and membranes such as the pleural and peritoneal layers, suggests a mechanism underlying the seed and soil hypothesis introduced by Stephen Paget in 1889 and yields realistic values for metastatic inefficiency. We propose a number of key experiments to test these concepts. The homeostatic pressure as introduced in this work could constitute a quantitative, experimentally accessible measure for the metastatic potential of early malignant growths.Comment: 13 pages, 11 figures, to be published in the HFSP Journa

Two-scale Moving Boundary Dynamics of Cancer Invasion:Heterotypic Cell Populations Evolution in Heterogeneous ECM

This book contains a collection of original research articles and review articles that describe novel mathematical modeling techniques and the application of those techniques to models of cell motility in a variety of contexts. The aim is to highlight some of the recent mathematical work geared at understanding the coordination of intracellular processes involved in the movement of cells. This collection will benefit researchers interested in cell motility as well graduate students taking a topics course in this area.

Elongation, proliferation & migration differentiate endothelial cell phenotypes and determine capillary sprouting

<p>Abstract</p> <p>Background</p> <p>Angiogenesis, the growth of capillaries from preexisting blood vessels, has been extensively studied experimentally over the past thirty years. Molecular insights from these studies have lead to therapies for cancer, macular degeneration and ischemia. In parallel, mathematical models of angiogenesis have helped characterize a broader view of capillary network formation and have suggested new directions for experimental pursuit. We developed a computational model that bridges the gap between these two perspectives, and addresses a remaining question in angiogenic sprouting: how do the processes of endothelial cell elongation, migration and proliferation contribute to vessel formation?</p> <p>Results</p> <p>We present a multiscale systems model that closely simulates the mechanisms underlying sprouting at the onset of angiogenesis. Designed by agent-based programming, the model uses logical rules to guide the behavior of individual endothelial cells and segments of cells. The activation, proliferation, and movement of these cells lead to capillary growth in three dimensions. By this means, a novel capillary network emerges out of combinatorially complex interactions of single cells. Rules and parameter ranges are based on literature data on endothelial cell behavior in vitro. The model is designed generally, and will subsequently be applied to represent species-specific, tissue-specific in vitro and in vivo conditions.</p> <p>Initial results predict tip cell activation, stalk cell development and sprout formation as a function of local vascular endothelial growth factor concentrations and the Delta-like 4 Notch ligand, as it might occur in a three-dimensional in vitro setting. Results demonstrate the differential effects of ligand concentrations, cell movement and proliferation on sprouting and directional persistence.</p> <p>Conclusion</p> <p>This systems biology model offers a paradigm closely related to biological phenomena and highlights previously unexplored interactions of cell elongation, migration and proliferation as a function of ligand concentration, giving insight into key cellular mechanisms driving angiogenesis.</p

Module-based multiscale simulation of angiogenesis in skeletal muscle

<p>Abstract</p> <p>Background</p> <p>Mathematical modeling of angiogenesis has been gaining momentum as a means to shed new light on the biological complexity underlying blood vessel growth. A variety of computational models have been developed, each focusing on different aspects of the angiogenesis process and occurring at different biological scales, ranging from the molecular to the tissue levels. Integration of models at different scales is a challenging and currently unsolved problem.</p> <p>Results</p> <p>We present an object-oriented module-based computational integration strategy to build a multiscale model of angiogenesis that links currently available models. As an example case, we use this approach to integrate modules representing microvascular blood flow, oxygen transport, vascular endothelial growth factor transport and endothelial cell behavior (sensing, migration and proliferation). Modeling methodologies in these modules include algebraic equations, partial differential equations and agent-based models with complex logical rules. We apply this integrated model to simulate exercise-induced angiogenesis in skeletal muscle. The simulation results compare capillary growth patterns between different exercise conditions for a single bout of exercise. Results demonstrate how the computational infrastructure can effectively integrate multiple modules by coordinating their connectivity and data exchange. Model parameterization offers simulation flexibility and a platform for performing sensitivity analysis.</p> <p>Conclusions</p> <p>This systems biology strategy can be applied to larger scale integration of computational models of angiogenesis in skeletal muscle, or other complex processes in other tissues under physiological and pathological conditions.</p

Ensemble Analysis of Angiogenic Growth in Three-Dimensional Microfluidic Cell Cultures

We demonstrate ensemble three-dimensional cell cultures and quantitative analysis of angiogenic growth from uniform endothelial monolayers. Our approach combines two key elements: a micro-fluidic assay that enables parallelized angiogenic growth instances subject to common extracellular conditions, and an automated image acquisition and processing scheme enabling high-throughput, unbiased quantification of angiogenic growth. Because of the increased throughput of the assay in comparison to existing three-dimensional morphogenic assays, statistical properties of angiogenic growth can be reliably estimated. We used the assay to evaluate the combined effects of vascular endothelial growth factor (VEGF) and the signaling lipid sphingoshine-1-phosphate (S1P). Our results show the importance of S1P in amplifying the angiogenic response in the presence of VEGF gradients. Furthermore, the application of S1P with VEGF gradients resulted in angiogenic sprouts with higher aspect ratio than S1P with background levels of VEGF, despite reduced total migratory activity. This implies a synergistic effect between the growth factors in promoting angiogenic activity. Finally, the variance in the computed angiogenic metrics (as measured by ensemble standard deviation) was found to increase linearly with the ensemble mean. This finding is consistent with stochastic agent-based mathematical models of angiogenesis that represent angiogenic growth as a series of independent stochastic cell-level decisions

3D Multi-Cell Simulation of Tumor Growth and Angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors

Adhesion Failures Determine the Pattern of Choroidal Neovascularization in the Eye: A Computer Simulation Study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression

Modeling Planarian Regeneration: A Primer for Reverse-Engineering the Worm

A mechanistic understanding of robust self-assembly and repair capabilities of complex systems would have enormous implications for basic evolutionary developmental biology as well as for transformative applications in regenerative biomedicine and the engineering of highly fault-tolerant cybernetic systems. Molecular biologists are working to identify the pathways underlying the remarkable regenerative abilities of model species that perfectly regenerate limbs, brains, and other complex body parts. However, a profound disconnect remains between the deluge of high-resolution genetic and protein data on pathways required for regeneration, and the desired spatial, algorithmic models that show how self-monitoring and growth control arise from the synthesis of cellular activities. This barrier to progress in the understanding of morphogenetic controls may be breached by powerful techniques from the computational sciences—using non-traditional modeling approaches to reverse-engineer systems such as planaria: flatworms with a complex bodyplan and nervous system that are able to regenerate any body part after traumatic injury. Currently, the involvement of experts from outside of molecular genetics is hampered by the specialist literature of molecular developmental biology: impactful collaborations across such different fields require that review literature be available that presents the key functional capabilities of important biological model systems while abstracting away from the often irrelevant and confusing details of specific genes and proteins. To facilitate modeling efforts by computer scientists, physicists, engineers, and mathematicians, we present a different kind of review of planarian regeneration. Focusing on the main patterning properties of this system, we review what is known about the signal exchanges that occur during regenerative repair in planaria and the cellular mechanisms that are thought to underlie them. By establishing an engineering-like style for reviews of the molecular developmental biology of biomedically important model systems, significant fresh insights and quantitative computational models will be developed by new collaborations between biology and the information sciences

Adrenergic regulation during acute hepatic infection with Entamoeba histolytica

Oxidative stress and transcriptional pathways of nuclear factor erythroid 2-related factor 2 (Nrf2) and nuclear factor kappa-B (NF-κB) are critically involved in the etiopathology of amebic liver abscess (ALA). In this work, we studied the relationship between the adrenergic nervous system and ALA in the hamster. ALA was visible at 12 h of infection. While 6-hydroxidopamine (6-OHDA) decreased infection, propranolol (β-adrenergic blocker) treatment was associated with less extensive liver damage, and phentolamine treatment (α-adrenergic blocker) significantly reduced ALA compared to 6-OHDA and propranolol. Serum enzymatic activities of alanine aminotransferase (ALT) and γ-glutamyl transpeptidase (γ-GTP) were increased at 12 h post-infection. Chemical denervation and α and β-adrenergic blockers decreased ALT to normal levels, while 6-OHDA and propranolol showed a trend to decrease γ-GTP but phentolamine significantly reduced γ-GTP. Amebic infection increased oxidized glutathione (GSSG) and decreased both reduced glutathione (GSH) and the GSH/GSSG ratio. Propranolol and 6-OHDA showed a tendency to decrease GSSG. However, GSH, GSSG and GSH/GSSG returned to normal levels with phentolamine. Furthermore, amebic infection increased pNF-κB and interleukin-1β (IL-1β), and showed a tendency to decrease hemoxigenase-1 (HO-1), but not Nrf2. Chemical denervation showed a trend to decrease pNF-κB and IL-1β, and neither Nrf2 nor HO-1 increased significantly. In addition, NF-κB and IL-1β were attenuated by propranolol and phentolamine treatments, although phentolamine showed significant overexpression of Nrf2 and HO-1. This suggests that the adrenergic system may be involved in oxidative stress and in modulation of the Nrf2 and NF-κB pathways during ALA development

Towards Predicting the Response of a Solid Tumour to Chemotherapy and Radiotherapy Treatments: Clinical Insights from a Computational Model

In this paper we use a hybrid multiscale mathematical model that incorporates both individual cell behaviour through the cell-cycle and the effects of the changing microenvironment through oxygen dynamics to study the multiple effects of radiation therapy. The oxygenation status of the cells is considered as one of the important prognostic markers for determining radiation therapy, as hypoxic cells are less radiosensitive. Another factor that critically affects radiation sensitivity is cell-cycle regulation. The effects of radiation therapy are included in the model using a modified linear quadratic model for the radiation damage, incorporating the effects of hypoxia and cell-cycle in determining the cell-cycle phase-specific radiosensitivity. Furthermore, after irradiation, an individual cell's cell-cycle dynamics are intrinsically modified through the activation of pathways responsible for repair mechanisms, often resulting in a delay/arrest in the cell-cycle. The model is then used to study various combinations of multiple doses of cell-cycle dependent chemotherapies and radiation therapy, as radiation may work better by the partial synchronisation of cells in the most radiosensitive phase of the cell-cycle. Moreover, using this multi-scale model, we investigate the optimum sequencing and scheduling of these multi-modality treatments, and the impact of internal and external heterogeneity on the spatio-temporal patterning of the distribution of tumour cells and their response to different treatment schedules